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ABSTRACT 

The compact dark objects with very large masses residing at the centres of galaxies 
are believed to be black holes. Due to the gravitational lensing effect, they would cast 
a shadow larger than their horizon size over the background, whose shape and size 
can be calculated. For the supermassive black hole candidate Sgr A*, this shadow 
spans an angular size of about 50 micro arc second, which is under the resolution 
attainable with the current astronomical instruments. Such a shadow image of Sgr A* 
will be observable at about 1 mm wavelength, considering the scatter broadening by 
the interstellar medium. By simulating the black hole shadow image of Sgr A* with 
the radiatively inefficient accretion flow model, we demonstrate that analyzing the 
properties of the visibility function can help us determine some parameters of the 
black hole configuration, which is instructive to the sub-millimeter VLBI observations 
of Sgr A* in the near future. 

Key words: black hole physics - relativity - methods: numerical - scattering - 
Galaxyicentre - sub-millimeter - techniquesrinterferometric. 



1 INTRODUCTION 

There is compelling evidence that a 4 x 10® black 
hole, associated with an extremely compact radio source 
Sagittarius A* (Ser A* ), resides at the Galactic Centre 



l|Schodel et al.ll2003l : IChez et al.l2005l : iMelia fc Falckel200ll ) . 

Recent Very Long Baseline Interferometry (VLBI) observa- 
tions reveal_th^tthB intrinsic size of Sgr A* is only about 1 
AU (|Shen et al.ll2005l ). For such a high degree of mass con- 
centration, the relaxation time scale is much shorter than 
the age of our Galaxy, and thus almost all the theoretical 
models predict the format ion of a sup ermassive black hole 
via gravitational collapse (|Maoj|l998l ). 

However, to conclusively prove that Sgr A* is a 
supermassive black hole requires observations close to its 
event horizon, where relativistic effects can not be ignored. 
According to General Relativity, the photon trajectory 
will be nd appreciably in a strong gravit at ional field (see, 
e.g., iMisner. Thorne fc Wheeled 1 19731 : IChandrasekhaJ 
1 19831 ). Because of this peculiar behavior, images of black 
holes and their immediate neighborhood we observe 
are highly distorted. The true physical structure in the 
very vicinity of the black hole must be decoded from 
the image via a ray t racing method (iFanton et al ] 1 19971: 
Falcke. MeUa fc Agoll |200(J: ISchnittman fc Bertschingerl 



20041 : ISchnittman. Krolik fc Hawlevll2006l ). In particular? if 



the impact parameter is smaller than a criti cal value, the 
phot on is doomed to fall into the black hole (iChandrasekhad 
Il983h . As seen by a distant observer, this effect allows the 
black hole to cast a shadow over the background source 
that is larger than its event horizon. 

The black hole candidate Sgr A*, at a distance of 8 kpc 
to us, spans an angular size of ~ 20 micro arc second (fias) in 
diameter in the sky, which is beyond the resolution attain- 
able with the current astronomical instruments. However, 
the shadow size is always of order WGM/c^ in diameter, 
corresponding to an angular size of 50^as for Sgr A*. In 
principle, this resolution can be achieved using VLBI tech- 
niques at a baseli ne length ^ 4 x 10^ km if obs erved at 1 mm 
wavelength (e.g., iFalcke. Meha fc Agolllioool ). 

In this paper, we adopt the ray-tracing method to map 
the Sgr A* black hole wi th the radiatively inefficient accre - 
tion flow (RIAF) model l|Yuan. Quataert fc Nara van"2003!). 
We then use our simulations to estimate the visibility func- 
tion to be observed by a VLBI array operating at, if can 
be designed in the future, sub-millimeter wavelengths. We 
are on the verge of resolving the shadow of the presumed 
supermassive black hole at the Galactic Centre. 

The paper is organized as follows. Section 2 briefly de- 
scribes the formalism we use for the ray-tracing algorithm. 
We show the simulation results in Section 3. Then we es- 
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timate the visibility functions in observations in Section 4. 
Finally we give a discussion in Section 5. 

Throughout this paper, we will work with geometric 
units where c = G = 1. 



2 HAMILTONIAN GEODESIC RAY-TRACING 
IN BLACK HOLE SPACETIME 

In a given spacetime metric, the trajectory of a photon 
is uniquely determined by the initial conditions. To effi- 
ciently simulate the image of a black hole and its environ- 
ment, we opt to impose our initial conditions at infinity, 
where the observer is, and integrate backwards in affine 
parameter toward the b lack hole (cf. iRauch fc BlandfordI 
[1994 ; Fanton ct all Il997l ). Starting from the standard ki- 
netic Lagrangian, we compute the su per- Ifamiltonian (cf. 
iMisner et all Il973l : iFeng fc Qiiil 120031 ) in Kerr metric but 
with spin parameter a = 0, since we discuss only the 
Schwarzschild black hole case in this work. Normal vari- 
ational approach yields the Hamiltonian equations of mo- 
tion in the phase space {q^,Pfj), where and p^j are 
Boyer-Linquist coordinate and the canonical momentum of 
the particle, respectively. If we choose the affine param- 
eter to be proper time per unit rest mass (which has a 
finite limit even when the rest mass is zero in the pho- 
ton case), then is interpreted as the covariant compo- 
nent of the four-momentum. Sim il ar pr oc ess can be seen in, 
e.g.. [ Schnittman fc Bertschingeij ()2004h : iBroderick fc Loebl 
(l200dV 

As discovered bv lBardeen et al.l (Il972l), an infalling pho- 
ton with its impact parameter ( Chandrasekharl Il983l ) less 
than \/27 M will never escape the vicinity of the black hole. 
As seen by an outside observer, the black hole will cast a 
shadow of radius rghadow = \/27M ~ 5.2M over the back- 
ground light source. As reported by many numerical cal- 
culations (e.g., Falcke et al, 2000; Takahashi 2004), rotating 
black holes will cast shadows of approximately the same size 
as well. 

In the following simulation, we specify initial conditions 
at a large distance (which we have chosen arbitrarily ro = 
200M), and then integrate backward toward the black hole 
to obtain the full geodesic. 



3 RESULTS FROM SIMULATION 

Fig. 1 is a sketch m ap, defining the inclination angle i 
l|Chandrasekharl 1 19831 ) and position angle Q used in our 
simulation. Three axes {X,Y,Z) point east, north and the 
direction aligning with the line of sight for an observer, re- 
spectively. The thick blue arrow represents the direction of 
the black hole's equatorial plane (the spin axis for the Kerr 
black hole). The inclination angle i £ [— 7r/2,7r/2] is the an- 
gle between the spin axis and the observer's line of sight 
(Z axis), and the position angle Q G [— 7r,7r] is the angle be- 
tween the projection of spin axis on X — Y plane and Y axis, 
positive by east and negative by west. We limit i G [0,7r/2] 
and 9 £ [~7r/2,7r/2] in the following discussion because of 
symmetry. 



3.1 Black Hole Shadow of Sgr A* with the 

Radiatively InefHcient Accretion Flow Model 

Even though there have been various theoretical models 
proposed t o explain the observations of Sgr A* (see, e.g., 
iMelia fc Fa lckc 2001), the precise structure of the dominant 
emitting region is still poorly known. In this paper, we onl y 
consider the RIAF model proposed bv lYuan et al.l l|2003l ). 
which can satisfy most observational results including the 
observed spectrum of Sgr A* from radio to X-ray, the flaring 
activity at both infrared and X-ray bands, and its extremely 
faint luminosity. 

In the RIAF model, the emission of Sgr A* from mil- 
limeter to sub-millimeter is dominated by the synchrotron 
radiation of thermal electrons. Thus one can determine the 
frequency-dependent emissivity j^, (r, 0) in the equatorial 
plane. Assuming a Gaussian distribution of electron density 
in the vertical direction of a RIAF, the emissivity distribu- 
tion function is then: 



j^(r,z) = >(r,0)exp(- — ) 



(1) 



where z is the vertical distance from the equatorial plane 
and H is the scale height determined by: 



H 



as 



(2) 



where as is the thermal sound speed and Qk is the Ke- 
plerian angular velocity. The absorption coefficient Kj^ is 
then given by Kirchhoff's Law. In this work, exactly the 
same RIAF model p a ramet ers in IVuan et al.l (|2003l ) and 
lYuan. Shen fc Huan3 (|2006l ) are adopted to preserve the 
good fit to the spectrum of Sgr A* and other observations. 

To construct an observed emission structure, we then 
perform a full radiation transfer calculation using the ray- 
tracing method outlined in Section 2. Here, the radiative 
transfer equation takes the form: 



_ . _ 

as 



(3) 



Since the frequency is red-shifted along the ge odesic and 
J./;/^ is a Lorentz invariant (see, e.g., IChan et al. 20061 : 
ISchnittman et al.ll2006l : lNoble et al]|2007l '). the observed spe- 
cific intensity is related to the emitted intensity by 



(4) 



where the gravitational red-shift factor is calculated to be 



5 = 



= ^ 1- 



2M 



(5) 



for a non-spinning black hole. We neglect the re-emission 
process when solving the radiation transfer equation. 

Furthermore, light rays emitted from the RIAF bend 
to us due to strong gravitation of the black hole, suffer ab- 
sorption and in addition, suffer interstellar scattering. That 
means, what we can observe is a scatter-broadened image. 
The two-dimensional scattering structure in the direction to 
the Galactic Centre is determined from fitting to the angu- 
lar size measurements as a function of the observing wave- 
length (cf. Shen et al.. ,2005;) . It is a Gaussian ellipse with 
full-widths at half maximum (FWHM) of major axis and 
minor axis in milli-arcsecond (mas) to be (1.39±0.02)A^ and 
(0.69± 0.06)A2 ( A measured in cm ), respectively and posi- 
tion angle ~ 80° l|Shen et al.ll2005l ). We take the conversion 
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that one gravitation radius corresponds to 5/ias angular size 
for Sgr A* with 4 x IO'^Mq mass at 8kpc distance. 

In Fig. 2, we show both the un-scattered (ray-tracing 
only) and scatter-broadened images at wavelengths of 1.3 
mm (first two columns of panels) and 3.5 mm (last two 
columns), respectively, obtained from our simulations with 
parameters = 0°, i = 7r/2,7r/4, (from top to bottom). 
For each two columns at two wavelengths, right one is ob- 
tained by convolving the left one (GR ray-tracing result) 
with the interstellar scattering. We use color to show the 
self-normalized specific intensity of the images. The abscissa 
axes point to west and the ordinal axes point to north. 
The non-scattered (intrinsic) images are different at difi'erent 
wavelengths. They exhibit predicted shadows as a result of 
relativistic effects. The scattered image becomes more ob- 
scured as the observing wavelength increases. It indicates 
that gradually with the availability of the high-resolution 
VLBI imaging at wavelengths of 1.3 mm or shorter, we will 
be at a very good position to unveil the shadow structure of 
the black hole. 



position angle (O), respectively. Thick solid lines are the 
upper and lower limits of the FMHMs with the best-fit 
scattering angular size of 1.39A^x0.69A^, and thick dashed 
lines are the upper and lower limits to account for the ilcr 
(0.02A^ X 0.06A^) in the scattering size. The measured an- 
gular sizes from high-resolution VLBA imaging (|Shen et al.l 
|2005i ) are indicated in Fig. 3 by the thin solid lines with 
the la u ncertainties by th e thin dotted lines. As already 
shown in I Yuan et al.l l|2006l ). within the uncertainties of the 
measurements and calculations, the predicted sizes from the 
RIAF model with i = 0° and = 0° are in reasonable agree- 
ment with the observations at two wavelengths. With a more 
thorough consideration shown in Fig. 3, it seems to suggest 
a geometry of the RIAF in Sgr A* to have a large inclina- 
tion angle i and a small position angle O. But it should be 
cautious since the results i n the mi nor axis are not accurate, 
especially at 3.5 mm fsee. lShen et al. 2005). It is clear to us 
that in order to get a reliable estimate of the geometry (i 
and 0) of the RIAF in Sgr A*, more accurate measurements 
are needed. 



3.2 Comparison with the VLBI Observations at 7 
and 3.5 mm 

Currently, VLBI observations can be steadily performed at 
millimeter wavelengths of 3.5 and 7 mm. Attempts to de- 
termine the Sgr A* structure with the VLBI observations, 
however, have suffered from the angular broadening caused 
by the diffractive scattering by the turbulent ionized inter- 
stellar medium, which dominates the resultant images with a 
A^-dependence apparent size. The development of the model 
fitting analysis by means of the amplitude closure relation 
along with the careful design of the observations at millime- 
ter wavelengths has greatly improved the accu racy of the 
size measurements of the observed image (e.g., IShen et al.l 
l2003l : iBower et aLll2004 IShen et al.ll2005l ). As a result, an 
intrinsic so urce size as comp act as 1 AU was first detected 
at 3.5 m m (IShen et al.l 20051 ). However, for the reason dis- 
cussed in lYuan et al.l (|2006l ). it is more meaningful to di- 
rectly compare the observed apparent VLBI image with the 
scatter-broadened image obtained from the simulations. The 
apparent images of Sgr A* can be depicted quantitatively by 
an elliptical Gaussian distribution with th e following param - 
eters of the major and minor-axis sizes l|Shen et 



4 VISIBILITY FUNCTIONS OF THE BLACK 
HOLE SHADOW IMAGES 



igparam- 
alT l2005h : 



0.724±0.001 mas by 0.384±0.013 mas and 0.2lt°;g^ mas by 
0.13lg j3 mas at wavelengths of 7 and 3.5 mm, respectively. 
The position angles at both wavelengths are about 80° , con- 
sistent wit h the orient a tion o f the scattering structure. 

While [Yuan et al] (|2006l ) only considered a special con- 
figuration, i.e., the RIAF is face-on {i = 0° and = 0°), here 
we will do a thorough investigation on the possible geome- 
try of the RIAF with respect to the presumed supermassive 
black hole in the Galactic centre based on the available VLBI 
measurements at both 3.5 and 7 mm. For this purpose, we 
first made a series of simulated after-scattering images with 
different geometric configuration (different combination of i 
and 0) between the black hole and the RIAF. Then, we fit 
the simulated images with an elliptical Gaussian distribu- 
tion to estimate the FWHMs of the major and minor axes. 

In Fig. 3, the FWHMs of the major and minor axes 
of the simulated scattering-broadened images of Sgr A* are 
plotted as a function of the inclination angle (i) and the 



As sh own in Fig. 2 (also see Fig. 1 in I Yuan. Shen fc Huangj 
|2006| ). the simulated apparent scatter-broadened images of 
Sgr A* at wavelengths of 1.3 mm and shorter cannot be well 
fitted by one elliptical Gaussian component, indicating that 
the interstellar medium scattering effects at 1.3 and 0.8 mm 
no longer dominate the observed morphology. So, it is very 
promising to image the shadow of the supermassive black 
hole Sgr A* with the high-resolution VLBI observations at 
sub- millimeter wavelengths. However, there is no such a sub- 
millimeter VLBI array available at present for the imaging 
observations. Therefore, we will not focus on the comparison 
of the images at wavelengths of 1 mm and shorter, but try 
the visibility analysis to show that some useful constraints 
on the emission of Sgr A* can be extracted from the limited 
visibility data without imaging. 

In the interferometric observation, the visibility func- 
tion V{u, v) measured on a baseline with coordinates (u, v) 
is related, by the Fourier transform, to the sky brightness 
distribution, i.e., the shadow image in our simulation. To 
provide instructive information to the real observation, we 
first perform two-dimensional Fourier transform to obtain 
the visibility function of shadows like those in Fig. 2. Then 
we analyze the visibility distribution along four specific di- 
rections in the sky plane with position angles of (i) 80° (E'), 
i.e., along the major axis of the scattering structure; (ii) 
— 10° (N'), i.e., along the minor axis of the scattering struc- 
ture; (in) 35° (NE'); and (iv) -55° (NW'). In Fig. 4, we 
show an example of such visibility slices with parameters, 
i = 78.75° and = -70° at A = 1.3mm. The abscissa is 
projected baseline length in the units of lO^A and the ordi- 
nate the normalized visibility amplitude. We plot slices in 
solid, long-dashed, short-dashed, and dotted line for E', N', 
NE', and NW' directions, respectively. To depict these visi- 
bility slices, we introduce two characteristic baseline lengths: 
E (marked with open circle) to denote the baseline length at 
which a normalized visibility decreases to 0.5, and a (marked 
with filled circle) to denote the baseline length at which the 
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visibility reaches its first valley. Such a minimum in the vis- 
ibility distribution is important because it implies that the 
image can no longer be described by a single Gaussian but 
must have an additional structure, which could be related to 
the black hole shadow. Visibilities at baselines longer than a 
have quite limited signal-to-noise ratio at millimeter wave- 
lengths. The typical detection limit with the current VLBA 
(Very Long Baseline Array) at 3.5 mm is about 0.12 mJy. 

It is obvious that different sets of geometrical param- 
eters (?" and &) will result in different visibility profiles at 
different wavelengths. For a given orientation of the RIAF 
model at a given wavelength, we can obtain two character- 
istic baseline lengths (E and a) along each of the four direc- 
tions (E', N', NE', and NW'), i.e. (Ee/, oe'), (Eat/, ctjv'), 
(EjvB', cTjvB')) S'lid (EivvK'i <^NW')- We further list these in 
the order of their lengths and denote them as E„ and On 
with n=l (the shortest) to 4 (the longest). 

Our simulations show that the scattering effect at 
3.5 mm is still quite large (see Fig. 2), thus we cannot get 
enough information of any fine structure of the shadow. To 
minimize the scattering effect, we consider the observations 
to be performed at other two shorter wavelengths 1.3 and 
0.8 mm. Shown in Fig. 5 arc the two characteristic baseline 
lengths (S and a) as a function of the inclination angle i 
at 1.3 (upper two panels) and 0.8 mm (lower two panels). 
With a specific i, each of S„ and a„ (n=l— 4) can vary 
within a region due to the different position angle Q. These 
allowed regions for Ei, E2, E3, and E4 (or, ai, 02, cr-i, and 
174) are shown as the darkest grey with long-dashed border 
lines, lightest grey with dotted border lines, second lighter 
grey with short-dashed border lines, and white with solid 
border lines, respectively. From these plots, we summarize 
some useful properties as follows: 

(i) Most of the characteristic baseline lengths arc in a 
range of (1—4) x 10^ km, which should be achievable un- 
der current VLBI observation conditions. Therefore, VLBI 
observations of Sgr A* at sub-millimeter wavelengths axe ex- 
pected to provide very important results in the near future. 
Some extreme cases may require a project baseline length 
longer than 5x10"^ km for (74. 

(ii) When the RIAF is nearly face-on [i as 0°), the four 
Es (or as) arc converged to a single similar number, suggest- 
ing that the four visibility profiles are almost identical and 
thus not orientation-dependent. This is mainly due to the 
isotropy of the intrinsic shadow image without scattering. 
The roughly East- West (80° in position angle) elongated 
scattering structure is very compact, ^ 15/ias along the ma- 
jor axis at about 1 mm wavelength, and thus will not affect 
the symmetry of the final scatter-broadened image. 

(iii) With a moderate inclination angle i, the differences 
between these projected baselines at different directions are 
quite significant. In fact, the four Es (or as) arc always range 
in four separate regions (as shown in Fig. 5) because of the 
asymmetry of the image. The possible region, caused by the 
different position angle G, for each E and a becomes larger 
when the inclination angle i increases. 

Therefore, in principle we would be able to constrain the 
inclination angle (i) of the RIAF in Sgr A* by comparing the 
above-mentioned characteristic baseline lengths determined 
from the future 1.3 and/or 0.8 mm VLBI experiments with 
the simulation demonstrated in this paper. It can be inferred 



from Fig. 5 that when the RIAF has a large inclination angle 
(in our simulation of Sgr A*, i >~ 30°), both Es and as 
become very sensitive to Q. 

No wonder, the position angle Q should have effects on 
the characteristic baseline lengths too. For clarity, wo define 
here another three normalized differences in the character- 
istic baseline lengths: 

S'l = (Ejv- -EbO/(S4-Ei), 

52 = (EjvB' — Ejviv')/(E4 — El), (6) 

53 = (Efi/ — Ejv£')/(^4 — El). 

Note that E4 and Ei represent the longest and shortest char- 
acteristic baseline lengths, respectively. If the RIAF is nearly 
face-on (i « 0°), the scatter-broadened image of Sgr A* 
is roughly E-W elongated, resulting in the E^y/ being the 
longest and E^/ the shortest (i.e., EAr/=E4 and E£;'=Ei), 
and EjvB'=Ejvw because of its symmetry. Therefore, it is 
always the case that Si = 1, 52 = and — 1 < 5,3 < for 
any Q. This means the face-on case is not sensitive to the 
position angle, consistent with the almost identical charac- 
teristic baseline lengths shown in Fig. 5 when i « 0° . 

Shown in Fig. 6 are plots of 5*1, S2 and S3 as a function 
of the position angle {&). At a given position angle 0, these 
values vary with the inclination angle i and have a range 
represented by grey region with dashed border lines (at 1.3 
mm) and dotted region with dotted border lines (at 0.8 mm). 

Assuming that the RIAF model is a physically reason- 
able description of the observed emission of Sgr A*, our 
simulation indicates that starting with the limited visibil- 
ity measurements on the four baselines in the directions 
mentioned above from the future sub-millimeter (1.3 and 
0.8 mm) VLBI experiments, we have a very good chance 
to get tight constraints on the configuration (i and G) of 
the RIAF around the central supermassive black hole of 
Sgr A*. That is, we can understand the geometry of the 
radio-emitting region surrounding Sgr A* in advance, even 
without imaging it with a completely (w, v) coverage of sub- 
millimeter/millimeter VLBI array. It is important to indi- 
cate that the detection of crs, i.e., valleys in the visibility 
profile (see Fig. 4 and Fig. 5) , are very important to prove 
that there is indeed a supermassive black hole in the centre 
of our Galaxy. 



5 DISCUSSION AND SUMMARY 

With the emissivity and absorption coefficient from the ra- 
diatively inefficient accretion flow model for the Galactic 
Centre black hole candidate Sgr A*, we solved the radia- 
tive transfer using the ray-tracing code to get the simulated 
shadow images of Sgr A* with an arbitrary set of geomet- 
rical parameters, i.e., i and G. By taking into account the 
interstellar scattering, we can obtain the apparent images at 
different wavelengths. 

Comparison of the simulated image sizes with those 
measurements at both 7 and 3.5 mm seems to prefer a RIAF 
with a large inclination angle and small position angle for 
Sgr A*. But, large uncertainties in both the size measure- 
ments and the interstellar scattering relationship make these 
very uncertain at the moment. Even though, our simulations 
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show that the future sub-millimeter VLBI images of Sgr A* 
would reveal a significant deviation from the single Gaussian 
structure, which can be used to probe the geometry of the 
RIAF of Sgr A*. 

Analysis of the scatter-broadened images also show that 
an observing wavelength of 1.3 mm or shorter is needed to 
identify the shadow shape. However, currently it is difficult 
to construct a VLBI array operated at 1.3 mm to obtain such 
an image. We therefore first perform the Fourier transform 
to obtain the corresponding visibility functions of images at 
1.3 and 0.8 mm, then analyze these visibility profiles along 
some specific directions. By introducing several characteris- 
tic baseline lengths (E„ and (t„) and the normalized differ- 
ences {S„), we demonstrate that it is possible to determine 
fundamental parameters (i and 6) of Sgr A* from the lim- 
ited detections on some baselines of VLBI observations at a 
wavelength of 1 mm or shorter (cf. Fig. 5 and Fig. 6). We 
can understand the geometry of the radio-emitting region 
surrounding Sgr A* even without imaging it with a com- 
pletely {u,v) coverage of a sub-millimeter/millimeter VLBI 
array. 

It should be mentioned that the four directions chosen 
in our analysis are somehow arbitrary. Any other combi- 
nations would also provide useful constraints on the struc- 
ture of Sgr A*. For example, it would be interesting to 
choose such two specific slices, one along the galactic plane 
and the other perpendicular to it. In practice, we can per- 
form the similar simulations with the projected baselines to 
be exactly the same as those from the real sub-millimeter 
VLBI experiments in the future. Of course the RIAF model 
adopted here is still quite simple. We s et the static blac k 
hole, use Pseudo-Newtonian potential (|Yuan et al.l |200^ . 
and simply assume a Gaussian distribution of the density as 
the vertical structure of the accretion fiow. And the possible 
dynamical role of an ordered magnetic field is not consid- 
ered. We also neglect the Doppler efi'ect of the radial and 
angular motions of the accretion flow in the radiative trans- 
fer calculation. However in principal, the visibility analysis 
outlined in this paper should be applicable, even if Sgr A* 
is a spinning black hole or with other accretion flow geome- 
tries and emission models. As a necessary flrst step trying 
to solve this complicated problem, this work definitely needs 
further efforts in both modelling and observations. And a 
more quantitative study will be required in future work. 

Furthermore, that the characteristic baseline lengths 
are about 10^ km and the expected observing wavelength 
is about 1 mm encourage us to propose the development 
of sub-millimeter VLBI array. The image of the vicinity of 
the predicted black hole in Sgr A* is expected to be more 
easily detected than any other super-massive black hole can- 
didates because of its closest distance to us. If the shadow 
image could be finally observed, it will be strong evidence 
for the General Relativity theory. 



Sciences, and sponsored by Program of Shanghai Subject 
Chief Scientist (06XD14024). Z.-Q. Shen and F. Yuan ac- 
knowledge the support by the One-Hundred- Talent Program 
of Chinese Academy of Sciences. 



REFERENCES 

Bardeen J.M., Press W.H., Teukolsky S.A., 1972, ApJ, 178, 
347 

Bower G. C, Falcke H., Herrnstein R. M., Zhao J., Goss 

W. M., Backer D. C. 2004, Sci, 304, 704 
Broderick A.E., Loeb A., 2006, MNRAS, 367, 905 
Chan C.K., Liu S., Fryer C.L., Psaltis D., Ozel F., Rocke- 
feller G., Melia F., 2006, astro-ph/0611269 
Chandrasekhar S., 1983, The Mathematical Theory of 

Black Holes, Oxford Univ. Press, Oxford 
Falcke H., Meha F., Agol E., 2000, ApJ, 528, 13 
Fanton C, Calvani M., Felice F., Cadez A. 1997, PASJ, 49, 
159 

Feng K., Qin M.-Z., 2003, Symplectic Geometric Algo- 
rithms For Hamiltonian System, Zhejiang Science & Tech- 
nology Press, Hangzhou. 

Ghez A.M., Salim S., Hornstein S.D., Tanner A. Lu J.R., 
Morris M., Beckhn E., Duchene G., 2005, ApJ, 620, 744 

Maoz E., 1998, ApJ. 494, L181 

Meha F., Falcke H., 2001, ARA&A, 39, 309 

Misner C.W., Thorne K.S., Wheeler J.A., 1932, Gravita- 
tion, W.H. Freeman & Company, San Francisco 

Narayan R., Yi I., 1994, ApJ, 428, L13 

Noble S.C., Leung P.K., Gammie C.F., Book L.G., 2007, 
|astro-ph/0701778 

Ranch K.P., Blandford R.D., 1994, ApJ, 421, 46 
Schnittman J.D., Bertschinger E., 2004, ApJ, 606, 1098 
Schnittman J.D., Krolik J.H., Hawley J.F., 2006, ApJ, 651, 
1031 

Schodel R., Ott T., Genzel R., Eckart A., Mouawad N., 

Alexander T., 2003, ApJ, 596, 1015 
Shakura N.I., Sunyaev R.A., 1973, A&A, 24, 337 
Shen Z.-Q., Liang M.C., Lo K. Y., Miyoshi M., 2003, As- 

tron. Nachr., 324, SI, 383 
Shen Z.-Q., Lo K.Y., Liang M.C., Ho P.T.P., Zhao J.H., 

2005, Nat, 438, 62 
Takahashi R., 2004, ApJ, 611, 996 

Yuan F., Quataert E., Narayan R., 2003, ApJ, 598, 301 
Yuan F., Shen Z.-Q., Huang L., 2006, ApJ, 642, L145 

This paper has been typeset from a TpjX/ ffl^jX file prepared 
by the author. 



ACKNOWLEDGMENTS 

We would like to thank two anonymous referees for their 
critical comments on the manuscript. This work was sup- 
ported in part by the National Natural Science Foundation 
of China (grants 10573029, 10625314, and 10633010) and the 
Knowledge Innovation Program of the Chinese Academy of 



6 L. Huang et al. 




Figure 1. Sketch map of the black hole. X-axis is to east, K-axis to north and Z-axis aligns with the line of sight of an observer. The 
blue arrow represents the direction of the black hole's equatorial plane which is aligned with the mid-plane of the RIAF (see text), and 
i and are the inclination angle and position angle, respectively. 
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Figure 2. Simulated images of tlie RIAF of Sgr A* witli © = 0° and i = it/2, it/A, (from top to bottom) at wavelengtlis of 1.3 mm (left 
two columns of panels) and 3.5 mm (right two columns of panels). For each two columns at two wavelengths, right one is obtained by 
convolving the left one (GR ray-tracing result) with the interstellar scattering. In each panel, the intensity is normalized to itself to show 
detailed structure. The units in the X- and Y-axes are gravitational radius Vg, which is 6x10 cm or bfias for Sgr A* with 4 X lO'^Afo 
mass at 8 kpc distance. 
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Figure 3. Comparison of the simulated sizes of Sgr A* with the real mcasmcments at 7 mm (upper four panels) and 3.5 mm (lower 
four panels). Panels in left column, from top to bottom, show sizes (i.e. FWHMs) of the scatter-broadened images along directions of 
major axis at 7 mm, minor axis at 7 mm, major axis at 3.5 mm, and minor axis at 3.5 mm, as a function of the inclination angle (i) of 
the RIAF, while the panels in right column show these sizes as a function of the position angle (0) of the RIAF. In each panel, the thick 
solid lines represent upper and lower limits of FWHMs of the images broadened by a constant scattering effect of 1.39A'^x0.69A'^ with a 
position angle of 80° . The upper and lower limits to the sizes are due to different geometrical structure, i.e. different position angles at a 
fixed inclination angle i, or different inclination angles at a fixed position angle 0. After considering the uncertainties in the scattering 
anerular sizes, the unner and lower limits are nlotted with the thick dashed lines. The measured sizes and errors are renresented bv those 
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Figure 4. The predicted visibility function of Sgr A* (in this case, at 1.3 mm, with i=78.75° and ©= —70°). The four curves represent 
the normalized visibility amplitude profiles at four position angles, i.e. 80° (solid line), —10° (long-dashed line), 35° (short-dashed line) 
and —55° (dotted line). We mark where a profile reaches 0.5 (open circle) and its first valley (filled circle). 
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Figure 5. The projected characteristic baseline Iciigtlis vcrsus^flic'Tncfination angle i at 1.3 mm (upper two panels) and 0.8 mm (lower 
two panels). With a specific inclination angle i, each of S„ and cr„ (n=l-4) can vary with the position angle 0. These allowed regions 
for El, E2, E3, and E4 (or, cti, erg, cs, and (T4) are shown as the darkest grey with long-dashed border lines, lightest grey with dotted 
border lines, second lighter grey with short-dashed border lines, and white with solid border lines, respectively. 
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Figure 6. The normalized difference of tlie cliaracteristic baseline lengths (Si, 52, and S3, from top to bottom) as a function of the 

position angle ©. For a given position angle (0), Si, S2, and S3 can vary with the inclination angle t, and have a range represented by 
the grey region with dashed border lines and the dotted region with dotted border lines for observations at 1.3 and 0.8 mm, respectively. 



